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CHAPTER  1 

EDGE-SOURCE  ACOUSTIC  GREEN’S  FUNCTION 
FOR  AN  AIRFOIL  OF  ARBITRARY  CHORD 


SUMMARY 

Approximations  are  derived  for  the  three-dimensional,  time-harmonic  acoustic  Green’s 
function  whose  normal  derivative  vanishes  on  the  surface  of  an  airfoil  of  finite  thickness  and 
chord  t  for  source  locations  in  the  neighborhood  of  either  the  leading  or  trailing  edge.  The 
acoustic  wavelength  is  assumed  to  be  large  relative  to  the  airfoil  thickness,  but  no  restriction 
is  placed  on  its  magnitude  relative  to  i.  A  multiple  scattering  calculation  is  performed 
for  high  frequencies  that  involves  the  successive  scattering  of  waves  from  the  leading  and 
trailing  edges  of  the  airfoil.  The  ‘principal  subseries’  of  the  expansion  is  summed  and  shown 
to  provide  an  excellent  approximation  for  the  Green’s  function  when  K0i  >  1,  where  k0  is 
the  acoustic  wavenumber.  The  solution  is  extended  down  to  k0£  =  0  by  interpolation  with 
the  corresponding  Green’s  function  for  an  airfoil  of  acoustically  compact  chord.  The  results 
extend  the  single  scattering  approximation  introduced  by  Amiet  [1],  and  are  illustrated 
by  application  to  the  problem  of  trailing  edge  noise  generated  by  nominally  steady,  low 
Mach  number  flow  past  the  airfoil.  Experiments  and  numerical  simulations  of  such  flows 
often  include  acoustic  frequencies  that  are  sufficiently  small  that  the  usual  assumption  of 
trailing-edge  noise  theory,  that  the  airfoil  is  semi-infinite,  is  not  valid. 
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1.1  INTRODUCTION 

The  ‘self  noise’  produced  by  nominally  steady,  high  Reynolds  number  flow  over  the 
trailing  edge  of  an  airfoil  is  usually  attributed  to  the  ‘diffraction’  of  turbulence  velocity 
fluctuations  by  the  edge  [2-13].  The  turbulence  is  generated  within  the  unstable  boundary 
layer  approaching  the  edge  and  is  swept  past  the  edge  by  the  mean  flow.  The  typical 
eddy  dimension  is  usually  small,  being  of  the  order  of  the  boundary  layer  thickness.  The 
sound  therefore  tends  to  be  of  relatively  high  frequency,  with  amplitude  and  spectral 
characteristics  dependent  on  both  the  geometry  of  the  airfoil  trailing  edge  and  on  the  details 
of  flow  separation  at  the  edge  [9  -  11,  14  -  16]. 

Most  estimates  of  trailing  edge  noise  have  involved  the  assumption  that  the  frequency 
is  sufficiently  large  that  the  acoustic  wavelength  may  be  regarded  as  very  much  smaller 
than  the  airfoil  chord.  Diffraction  theory  predictions  were  therefore  made  by  replacing  the 
airfoil  by  a  semi-infinite,  rigid  plane  [2  -  6].  Adamczyk  [17],  Amiet  [18,  19]  and  Martinez  & 
Widnall  [20]  have  studied  theoretically  the  sound  produced  by  turbulence  and  discrete  gusts 
incident  on  a  thin  plate  airfoil,  including  a  ‘single  scattering’  correction  to  allow  for  the 
influence  of  finite  chord;  the  unsteady  lift  was  first  calculated  and  the  radiation  subsequently 
determined  by  means  of  Curie’s  [21]  surface  integral  representation  of  aerodynamic  sound. 
A  similar  procedure  [22,  23]  has  been  proposed  for  including  the  influence  of  finite  chord  on 
trailing  edge  noise.  Such  corrections  are  usually  believed  to  be  small  at  the  high  frequencies 
spanned  by  the  edge  noise  spectrum.  However,  this  may  not  be  universally  true  at  the 
very  low  Mach  numbers  encountered  in  underwater  applications.  In  such  cases  it  is  often 
required  to  be  able  to  predict  the  edge  noise  from  data  derived  from  turbulence  wall  pressure 
measurements  near  the  edge  or  from  numerical  simulation  of  the  edge  flow  [15,  16].  This 
data  will  essentially  describe  the  hydrodynamic  (incompressible)  properties  of  the  motion 
near  the  edge,  and  cannot  be  used  to  predict  the  radiation  using  Curie’s  formula,  because 
the  latter  requires  accurate  phase  information  of  the  surface  pressure  in  the  acoustic  domain 
when  the  surface  is  not  acoustically  compact  [15]. 

This  conclusion  has  important  implications  for  the  validation  of  trailing  edge  noise 
predictions  obtained  by  numerical  simulations  at  low  Mach  numbers.  Current  numerical 
algorithms  are  restricted  to  relatively  low  Reynolds  numbers,  for  which  the  edge  noise 
spectrum  extends  down  to  frequencies  where  the  acoustic  wavelength  is  comparable  to  the 
airfoil  chord.  Similarly,  the  airfoil  thickness  near  the  trailing  edge  is  often  too  large  to 
permit  the  acoustic  radiation  to  be  estimated  with  confidence  from  thin  airfoil  theory  [8  - 11, 
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24  -  26].  Separation  near  the  edge  can  lead  to  the  production  of  large  scale  flow  structures 
that  generate  significant  acoustic  energy  at  low  frequencies.  In  these  circumstances  a  proper 
comparison  with  experiment  of  predictions  of  the  sound  from  the  numerical  simulations  can 
be  made  only  when  both  the  finite  chord  and  finite  thickness  of  the  airfoil  are  properly 
incorporated  into  the  aeroacoustic  theory.  That  is  the  subject  of  this  chapter. 

The  radiated  sound  can  be  calculated  in  terms  of  the  hydrodynamic  edge  flow  and  a 
suitable  Green’s  function  G  [14].  In  the  following  attention  is  confined  to  a  rigid  airfoil,  for 
which  Green’s  function  is  required  to  have  vanishing  normal  derivative  on  the  airfoil  surface. 
The  mean  flow  Mach  number  is  assumed  to  be  very  small  (appropriate  to  underwater 
applications)  so  that  the  acoustic  wavelength  may  always  be  regarded  as  large  compared 
to  the  airfoil  thickness.  An  approximate  representation  of  G  will  be  derived  in  two  steps. 
First  the  airfoil  chord  is  assumed  to  exceed  the  characteristic  wavelength  of  the  sound,  and 
the  corresponding  approximation  for  G  is  obtained  by  considering  a  multiple-scattering 
problem  in  which  sound  waves  are  successively  scattered  from  the  leading  and  trailing  edges 
of  the  airfoil.  This  generalizes  Amiet’s  [1]  single  scattering  approximation,  and  incorporates, 
in  addition,  the  influence  of  the  finite  thickness  of  the  edge;  it  is  also  more  convenient 
than  Amiet’s  original  formulation  because  it  does  not  require  an  explicit  evaluation  of  the 
unsteady  lift  distribution,  nor  does  it  require  the  aeroacoustic  (turbulence)  sources  to  be 
convecting  in  a  prescribed  manner.  Second,  this  high  frequency  approximation  is  used  in 
conjunction  with  Green’s  function  for  an  airfoil  with  an  acoustically  compact  chord  to  derive 
an  interpolated,  composite  Green’s  function  that  is  expected  to  be  valid  for  all  frequencies. 

The  relevant  integral  representations  of  trailing  edge  noise  at  low  Mach  numbers  are 
recalled  in  §2,  and  the  generalized  Green’s  function  is  derived  in  §3.  Application  is  made  in 
§4  to  determine  the  influence  of  finite  chord  on  traditional  predictions  of  trailing  edge  noise 
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1.2.  REPRESENTATION  OF  THE  EDGE  NOISE 

Consider  the  sound  generated  by  low  Mach  number,  turbulent  flow  of  a  homogeneous 
fluid  of  mean  density  p0  and  sound  speed  c0  in  the  vicinity  of  the  trailing  edge  of  a 
rigid,  two-dimensional  airfoil.  The  discussion  will  be  framed  in  terms  of  the  canonical, 
parallel-sided  airfoil  shape  illustrated  schematically  in  Figure  1,  which  has  chord  t  and 
uniform  thickness  h,  and  a  rounded  nose.  Blake  and  his  coworkers  [8-11]  have  made 
extensive  measurements  of  the  edge  flow  and  sound  produced  by  airfoils  of  this  type  with  a 
variety  of  rounded,  ‘knuckled’  and  beveled  trailing  edges.  Low  Mach  number  turbulent  flow 
past  this  airfoil  type  is  also  the  subject  of  numerical  simulations  by  Wang  [24  -  26]. 

Introduce  the  rectangular  coordinate  system  (xi,xi,x3)  with  origin  O  at  a  convenient 
point  near  the  trailing  edge  such  that  the  ‘upper’  and  ‘lower’  planar  surfaces  of  the  airfoil 
are  at  x%  —  ±|h,  the  rzq-axis  is  in  the  direction  of  the  mean  flow,  and  x$  is  parallel  to  the 
airfoil  span  (out  of  the  plane  of  the  paper  in  Figure  la).  The  mean  flow  velocity  outside  the 
boundary  layers  on  the  airfoil  is  assumed  to  be  at  sufficiently  low  subsonic  speed  U  (in  the 
positive  aq-direction)  that  the  convection  of  sound  may  be  neglected.  This  will  be  the  case 
provided  the  Mach  number  M  =  U/c0  <C  1;  this  condition  also  ensures  that  p0  and  c0  may 
be  regarded  as  constant  throughout  the  flow. 

When  these  conditions  are  satisfied  the  far  field  acoustic  pressure  fluctuations  p(x,  u))e~tut 
of  frequency  to  produced  by  the  interaction  of  the  turbulence  with  the  airfoil  can  be 
expressed  in  the  form  [14,  15] 

'  (fi  Av)(y,o;)d3y  -  v  ^  £2(y,u;)  A  ^(x,y,u;)  •  n  dS(y),  (1) 

where  v  is  the  fluid  velocity,  £2  =  curl  v  is  the  vorticity,  and  u  is  the  kinematic  coefficient  of 
shear  viscosity.  The  second  integral  is  taken  over  the  surface  S  of  the  airfoil  (with  surface 
element  ndS(y),  and  with  the  unit  normal  n  directed  into  the  fluid),  and  represents  the 
contribution  to  the  radiation  from  the  unsteady  skin  friction;  it  can  usually  be  ignored 
when  the  Reynolds  number  is  large.  The  function  G(x,  y,  u>)  is  the  time  harmonic  Green’s 
function.  This  has  outgoing  wave-behavior,  vanishing  normal  derivative  dG(x,  y;  cj) /dyn  =  0 
on  S,  and  satisfies 

(V2  +  K;)G  =  J(x-y),  (2) 

where  k0  =  oj/c0  is  the  acoustic  wavenumber. 

Equation  (1)  is  strictly  valid  when  the  flow  is  homentropic  [27],  when  the  pressure  p  can 
be  regarded  as  a  function  of  the  density  p.  It  can  be  used  to  calculate  the  sound  radiated 
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from  unsteady  (vortical)  regions  provided  the  source  terms  (ft  A  v)  (y,  w)  and  fi(y,w) 
are  known,  either  from  measurement  or  from  a  preliminary  numerical  simulation  of  the 
flow.  The  notation  here  implies  that  a  source  term  F(y,t),  say,  is  known  as  a  function  of 
source  position  y  and  of  the  time  t,  and  that  F( y,  w)  =  ^  F( y,  t)elujtdt  is  the  Fourier 

transform  with  respect  to  the  time. 

The  contribution  from  the  volume  integral  in  equation  (1)  can  be  recast  at  very  low  Mach 
numbers  as  a  surface  integral  over  S  involving  the  ‘upwash’  velocity  induced  on  S  by  the 
turbulence  [15].  This  representation  for  the  particular  case  of  an  airfoil  of  zero  thickness 
(, h  -  0)  corresponds  to  the  original  diffraction  theory  of  trailing  edge  noise  developed  by 
Chase  [3,  5]  and  Chandiramani  [4],  which  can  also  be  expressed  in  the  form 

p(x,  u>)  =  -  f  G(x,  y,  u>) Vpi(yi,  0,  y3, u)  ■  ndS(y). 

J  s 

=  ~  Loo  dy 3  l-i  °’  y3)  ^  tG^X’  y>  dyi  ^ 

where 

[G(x,  y,  w)]  =  G(x,  pi,  +0,  y3,  u)  -  G(x,  yu  -0,  y3,  u)  (4) 

is  the  jump  in  the  value  of  G  across  the  airfoil.  The  pressure  pi  is  equal  to  that  which  would 
be  produced  by  the  turbulence  if  the  airfoil  were  imagined  to  be  temporarily  removed;  it 
can  be  expressed  in  terms  of  the  boundary  layer  ‘blocked  surface  pressure’,  and  equation 
(3)  represents  the  edge  noise  in  terms  of  the  scattering  of  this  pressure  by  the  edge. 
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1.3.  GREEN’S  FUNCTION 
1.3.1  Zeroth  order  approximations 

At  low  Mach  numbers  the  trailing  edge  interaction  noise  tends  to  be  dominated  by 
contributions  to  the  integrals  in  (1)  from  the  immediate  neighborhood  of  the  edge  region. 
The  wavelength  of  the  sound  ~  5/M,  where  5  is  a  length  that  characterizes  the  eddy  size 
near  the  edge,  and  may  in  practice  be  of  the  order  of  the  boundary  layer  thickness  or  the 
thickness  h  of  the  airfoil.  In  either  case  the  acoustic  wavelength  greatly  exceeds  h  when 
M<1,  which  implies  that  the  turbulence  responsible  for  the  edge  noise  is  always  very 
much  closer  than  an  acoustic  wavelength  from  the  edge.  Then  the  nondimensional  source 
distance  K0\Jy \  +  y\  ~  \jy\  +  y\ / acoustic  wavelength  <  1.  When  the  observation  point  x 
is  in  the  acoustic  far  field  and  y  is  within  the  source  region,  Green’s  function  can  then  be 
expanded  in  terms  of  this  small  parameter: 

G(x,  y,  u)  =  G0(x,  y,  u)  +  Gi  (x,  y,  u)  +  . . . .  (5) 


In  the  particular  case  of  an  airfoil  whose  chord  £  is  acoustically  compact  (kq£  <C  1)  we 
find  for  K„yfyi  +  y\  -C  1  [14] 


Go(x,y,w) 


_et/c0|x-j/3i3|  iKoSm'tpsmO  y/iip*(y)elKo^  J'3l3l 

71 - — ,  Gi(x,y,w)  = - - 

47T  x  —  y3i3  47r  x-j/3i3| 


(6) 


|x  -  2/313!  -»>  00, 

where  i3  is  a  unit  vector  parallel  to  the  airfoil  span  (in  the  positive  ^-direction).  If 
r  =  \Jx\  +  x\  denotes  the  observer  distance  from  the  edge  of  the  airfoil  (the  a^-axis),  then 
(r,  9)  are  the  polar  coordinates  of  the  observer  position  relative  to  a  plane  rr3  =  constant, 
i.e.  (£1,2:2)  —  r(cos0,sin0);  =  sin-1(r/|x|)  is  the  angle  between  the  observer  direction  x 
and  the  edge  (see  Figure  lb).  The  function  (p*(y)  =  <p*(yi, 2/2)  of  the  source  position  y  is  a 
solution  of  Laplace’s  equation  that  may  be  interpreted  physically  as  the  velocity  potential 
of  ideal,  incompressible  flow  around  the  edge  (in  the  anticlockwise  direction),  and  therefore 
depends  on  the  geometrical  shape  of  the  edge.  At  source  distances  \]yl  +  y"i  from  the  edge 
that  are  large  compared  to  the  airfoil  thickness  h  but  small  compared  to  the  chord  £,  </?*(y) 
must  tend  to  the  following  expression  for  the  potential  of  flow  around  a  rigid  half -plane: 

<p* (y)  -*  \fr’ sin(#' /2),  h<^r'  <^£,  where  (2/1, 2/2)  =  r' (cos 9', sin 6').  (7) 
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In  the  opposite  extreme  in  which  k01  —>■  oo,  when  the  airfoil  chord  is  semi-infinite,  the 
zeroth  order  component  Go(x,  y,  to)  in  the  expansion  (5)  is  identical  with  that  in  (6)  for  the 
compact  airfoil,  but  [14,  15] 


<?i(x,y,w)  = 


— y/^sin^  ^sin  </?*(y)eIK°lx  2/3l3l  , - 

/  ..  ",  i  j  t  OO,  Ho'J V\  2/2  ^ ‘ 

7rv27rz|x  —  2/3I3 1 


(8) 


In  both  limits  Go  represents  the  sound  produced  by  a  point  source  at  y  when  scattering  by 
the  airfoil  is  neglected.  The  component  G\  gives  the  first  correction  due  to  the  presence  of 
the  edge  in  the  nearfield  of  the  source,  and  supplies  the  leading  approximation  to  the  edge 
noise  when  used  in  equations  (1)  or  (3). 


Without  loss  of  generality  it  will  be  assumed  in  what  follows  that  the  turbulence  edge 
sources  are  located  in  the  neighborhood  of  the  origin  0.  The  notation  may  then  be  simplified 
at  the  outset  by  setting  y3  =  0  in  the  equations  (6)  and  (8).  It  may  also  be  noted  that 
sin  ip  sin  6  =  cos  0  in  (6),  where  ©  is  the  angle  between  the  radiation  direction  x  and  the 
nominal  airfoil  normal  (the  positive  a^-axis).  In  the  compact  limit  the  term  G\  accordingly 
represents  an  acoustic  dipole  orientated  in  the  direction  of  the  mean  lift,  in  accordance  with 
Curie’s  [21]  theory. 


1.3.2  Single  scattering  approximation 

To  bridge  the  extensive  frequency  interval  between  the  approximations  (6)  and  (8) 
occupied  by  airfoils  of  large  but  finite  chord  l  we  introduce  the  notation 

G{x.,y,co)  =  Gi(x,y,w)  +  GleC^Y^)  +  Gte(x)Y,  w).  (9) 

Here  and  henceforth  the  zeroth  order,  monopole  component  Go  is  discarded,  since  it  plays 
no  role  in  the  production  of  edge  noise.  The  components  Gle>  Gte  will  respectively 
represent  additional  contributions  to  G  that  may  be  associated  with  scattering  at  the 
leading  and  trailing  edges  of  the  airfoil.  These  terms  arise  in  the  following  way. 

For  an  airfoil  of  large  chord  we  take  G\  in  equation  (9)  to  be  defined  by  equation  (8) 
(with  2/3  =  0)  for  an  airfoil  with  semi-infinite  chord.  This  represents  an  acoustic  disturbance 
whose  wavelength  is  much  larger  than  the  airfoil  thickness  {n0h  <  1)  propagating  as  a 
function  of  x  in  the  presence  of  an  airfoil  of  semi-infinite  chord.  It  therefore  produces 
scattered  waves  at  the  leading  edge,  which  are  in  turn  scattered  at  the  trailing  edge  and 
subsequently  rescattered  at  the  leading  edge,  etc.  The  functions  GLe  and  Gte  are  defined 
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to  account  for  the  aggregate  contributions  from  the  leading  and  trailing  edges  of  all  of  these 
scattering  events. 

Because  K0h  <C  1,  multiple  scattering  at  the  edges  can  be  calculated  by  replacing  the 
airfoil  by  one  of  zero  thickness.  Furthermore,  as  implied  in  the  previous  paragraph,  the 
scattering  calculation  for,  say,  the  leading  edge  is  performed  by  assuming  the  zero-thickness 
airfoil  to  extend  infinitely  far  downstream  from  the  edge  (to  be  semi-infinite );  the  resulting 
scattered  field  is  then  corrected  at  the  trailing  edge  by  a  similar  calculation  performed, 
again,  for  a  semi-infinite  airfoil.  Repeated  application  of  this  procedure  and  summation  over 
all  scattered  contributions  from  the  leading  and  trailing  edges  will  then  yield  GL E,  Gte- 
The  method  is  identical  to  that  used  by  Schwarzschild  [28]  for  the  simpler,  two-dimensional 
problem  of  the  diffraction  by  a  slit  of  a  normally  incident  plane  wave,  and  is  also  discussed 
by  Landahl  in  connection  with  two-dimensional  gust  loading  of  an  airfoil  [29]. 

The  calculation  can  be  carried  out  in  full  in  terms  of  an  expansion  parameter  ~  1/ k0£  simp, 
which  is  required  to  be  sufficiently  small.  The  resulting  multiple  scattering  expansion 
turns  out  to  be  analytic  in  k0  in  the  region  | /c0 1 ^  sin  *0  >  0.125  of  the  upper  complex  plane, 
which  suggests  that  the  initial  high  frequency  approximation  can  be  analytically  continued 
to  more  modest  frequencies.  On  this  basis  we  shall  postulate  the  existence  of  a  smooth 
interpolation  with  the  ( compact )  very  low  frequency  approximation  of  equation  (6).  The 
details  of  the  calculation  are  straightforward,  and  will  be  illustrated  by  consideration  of  the 
leading  edge  scattered  field  Gle- 


The  zeroth  order  approximation  G\  (equation  (8)  with  y 3  —  0)  is  first  expressed  as  the 
following  Fourier  integral  on  the  upper  and  lower  surfaces  of  the  zero  thickness  airfoil  in  the 


vicinity  of  the  leading  edge  at  Xi  = 


Gi(x,  y,uj)  = 


T^sina  ip  sin  (§)  <£*(y)em°|x| 


*v®|x| 


— sgn(x2)y?*(y)e14 
(27r)2\/7r 


dkdkz, 

V k0  +  k 


R o  =  \[k?0  -  Kot  >  1,  (10) 


where  sgn^)  =  ±1  according  as  rr2  =  ±0,  and  the  radiation  condition  requires  that  the 
imaginary  part  of  R0  =  sJk2q  -  h%  should  be  positive. 

The  interaction  of  G\  with  the  leading  edge  of  the  airfoil  may  now  be  cast  in  terms  of 
a  diffraction  problem  for  each  Fourier  component  proportional  to  sgn(:r2)et(fcxi+fc3X3).  The 
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diffracted  field  can  be  determined  exactly  when  the  airfoil  is  regarded  as  semi-infinite  (when 
the  trailing  edge  is  moved  to  xi  =  +00).  The  routine  details  of  the  calculation  are  discussed 
at  great  length  by  Noble  [30],  and  it  is  sufficient  here  to  quote  the  final  form  of  the  scattered 
field.  Denoting  this  first  scattered  contribution  from  the  leading  edge  by  G[ E,  we  find 


G[ E(x,  y,  u) 


— sgn(x2)^*(y)e  '* 
(2tt)20F 


® j ^3)  \&2 1  k,ty 

y/R0  +  k\jR0  —  ki{ki  —  k  +  iO) 


dkidk3dk, 


(11) 


where  x[  =  X\  +  £,  and  the  branch  of  7(^1,  k$)  —  y/c2  —  k\  —  k$  with  positive  imaginary 
part  is  taken  to  satisfy  the  radiation  condition. 


The  next  step  is  to  consider  the  interaction  of  this  wave  with  the  trailing  edge.  If, 
however,  all  such  higher  order  interactions  are  neglected  we  shall  obtain  a  single  scattering 
approximation  to  G(x,  y,ca).  This  corresponds  to  the  approximation  used  by  Amiet  [1,  18, 
19],  and  we  shall  pause  to  obtain  its  explicit  representation  when  the  observer  position  x 
moves  to  the  acoustic  far  field  of  the  airfoil.  This  is  done  by  first  evaluating  the  integrations 
with  respect  to  ki  and  k$  in  (11)  by  the  method  of  stationary  phase  [14],  which  supplies 


G[E(x>y^) 


y/2K^ip*{y)  sin  2  ip'  cos  (y)  ei(Ko|x'|+ **) 


x 


/: 


(27r)2y/7r|x,| 

a/aCq  sin  ip'  —  k  e~Mdk 


ac0|x'|  — >  00, 


(12) 


-00  \/7i0  sin  ip'  +  k(n0  sin  ip'  cos  9'  —  k  +  iO)  ’ 
where  x'  =  (xi  +  x2,  £3)  and  the  corresponding  angles  ip',  6'  are  defined  by  analogy  with 

the  angles  ip,  6  of  Figure  lb,  but  with  the  origin  at  the  leading  edge. 


This  formula  is  simplified  further  by  noting  that,  when  |x|  l  we  can  replace  x'  by 
x  in  the  denominator,  and  the  angles  ip'  ip,  9'  -¥  9.  When  k0£  is  large  the  principal 
contribution  to  the  remaining  integral  is  from  the  neighborhood  of  the  branch  point  of  the 
integrand  at  k  =  —  k0  sin  ip'  — *  —  K0sin^.  We  can  therefore  set  'ipi~^k  =  y/2n0  sin  ip 

in  the  numerator,  but,  apart  from  replacing  ip',  9'  by  ip,  9,  the  denominator  must  be 
left  intact,  because  both  terms  can  be  arbitrarily  small  near  the  branch  cut  when  9  — »  ir. 
In  this  limit  the  pole  of  the  integrand  produces  the  forward  scattered  wave  needed  to 
cancel  the  incident  wave  (10)  on  the  upstream  extension  of  the  airfoil.  However,  it  is  not 
necessary  to  make  any  further  approximations,  because  when  sin^'  —  k  is  replaced 
by  y/2 /c0  sin  ip  the  integral  can  be  evaluated  in  closed  form  in  terms  of  the  error  function 
«fw  =  *  Jq  e  When  the  result  is  combined  with  the  zeroth  order  term  (10)  we 

obtain  the  following  single  scattering  approximation  to  the  Green’s  function 
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G(x,y,u)  ~  Gi(x,y,o;)  +  G[E(x,y,o;) 


TTV^ih 


_  gtKo(|x#|— €sin^cos^) 


1  -  erf  e  ^  \  2k0£  sin  0  cos2 


I ,  /C0|x|  — »■  oo. 

(13) 


For  purposes  of  computation  the  error  function  of  complex  argument  may  be  expressed  in 
terms  of  the  real- valued  Fresnel  integral  auxiliary  functions  f(x),  g(x)  defined  in  §7.3  of 
[31],  in  which  case  the  single  scattering  approximation  becomes 


G(x,y,w) 


[sin 

7r\/27ri|x|  _  2 


7T 


Ko  X 


oo. (14) 


where  T(x)  =  g(x)  +  if(x). 


1.3.3  Multiple  scattering 

Return  now  to  the  multiple  scattering  calculation  of  Gle-  First  order  scattering  at  the 
leading  edge  produces  the  acoustic  field  G\ /E  given  by  (11).  A  diffraction  calculation  must 
be  performed  to  determine  the  field  produced  when  this  wave  is  scattered  back  to  the 
leading  edge  from  the  trailing  edge.  To  do  this  the  representation  (11)  is  first  used  to 
simplify  the  expression  for  the  first  order  wave  incident  on  the  trailing  edge  by  replacing 
the  integrations  with  respect  to  k  and  ki  by  the  corresponding  leading  order  term  in  the 
asymptotic  expansion  for  k,0£  1,  which  is  readily  found  in  the  form 


GxE(x,y,w) 


-sgn(z2V*(y)e- 

ttV2^ 


ITT 

-  [°°  y/R~0  c  e^'+^dk 3,  x2  =  ±0,  (15) 

J -OO 


where 


J2iK0£ 


t  = 


(16) 


2tt  iR0£ 

A  diffraction  calculation  is  then  performed  for  each  Fourier  component  proportional  to 
ei(Roxi+k3x3)  incident  on  the  trailing  edge,  as  before.  The  net  result  of  these  calculations  is 
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just  the  first  term  C?TE(x,y,u>),  say,  in  the  sequence  of  waves  scattered  from  the  trailing 
edge.  Just  above  and  below  the  airfoil  this  wave  is  found  to  be  given  by 


GxE(x,y,w)  = 


— sgn(x2)y?*(y)e 
(27r)2v/7r 


17T 

-SSI 


2&0e  e^Xl_*“k3X3) 


oo  \J k0  +  k(k  —  k o  —  zO) 


dkdks,  X2  =  ±0.  (17) 


The  integrand  can  be  simplified  by  replacing  k  —  R0  —  iO  by  —2 R0  when  X\  ~  —  £,  because 
when  K0i  is  large  the  main  contribution  from  the  integration  with  respect  to  k  is  from  the 
vicinity  of  the  branch  point  at  k  —  -R0.  When  the  result  is  combined  with  (10)  we  then 
find,  near  the  leading  edge, 


Gi(x,y,w)  +  G^E(x,y,w) 


-sgn(x2)9*(y)eI4  j-j 'o° 

(27r)2V7r  J  J — 00  \Ac<>  +  k 


[1  —  e]  dkdkz . 


(18) 


It  is  now  clear  that  repeated  application  of  the  calculations  leading  to  this  result  will 
produce  the  following  expansion  for  the  whole  of  the  field  incident  on  the  leading  edge 


Gi(x,y,w)  +  GxE(x,y,w)  +  G|E(x,y,a;)+  ••• 


Gi(x,y,w)  +  GTE(x,y,w) 


-sgn(a:2)¥>*(y)e,4  ^ooe^>+M 

(2tt)2v^  77-00  VR^Tk  nV  }  dMh' 


X\  ~  —£,  X2  =  ±0.  (19) 


Now  e  =  e2iR°e /2iriR0£,  so  that  the  series  converges  absolutely  for  \R0\£  >  1/27T  for  real 
values  of  R0,  and  everywhere  in  Im  R0  >  0  outside  a  region  enclosing  the  origin  bounded  by 
the  curve  2ir\R0\£  =  exp  (— 21m  R0£).  In  the  region  of  convergence 


OO 


OO 


E  M“  =  E 


-e 


2iR,0£  \  n 


2mK0£ 


e2iK0i  \  1 

2mR0£  J  ’ 


(20) 


so  that  equation  (19)  becomes 
G\  (x,  y,  u)  +  GTE(x,  y,  u)  « 


-sgn(o;2)y>*(y)e,4 

(27r)2v/7r 


If. 


gi(kxi+ksxz) 


00  a/aCq  +  k(  1  +  e2m°e /2-KiK0£) 


dkdk 


(21) 


The  component  GEe  of  Green’s  function  is  produced  by  diffraction  of  this  field  at  the 
leading  edge.  Performing  the  diffraction  calculation  in  the  manner  already  discussed,  and 
considering  the  limiting  form  in  the  acoustic  far  field,  where  «0|x|  1,  we  accordingly  find 
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Gle^Yiw) 


?r\/27ri|x|(l  +  e2iK^sin^ /2mK0£simp) 


1  —  erf  e  "  \  2k01  sin  ip  cos2 


e 


v^sin^y*(y)eiK°(lx,l^sin^  / 
mt$|x|(1  +  e2iKotsin'1’ /2niK0£ sin  ip)  V 

(22) 


I  k0£  sin  ip  cos2  | 


7T 


«0|x|  — »•  00. 


Similarly,  the  trailing  edge  component  GTE(x,y,w)  is  produced  by  the  diffraction  of 
G'le(x,  y,  w)  at  the  trailing  edge,  and  is  found  by  the  same  procedure  to  be  given  in  the 
acoustic  far  field  by 


GTE(x,y,w) 


_  (p*  (y  )eiKo{|x|+^  sin  V>(l+cos  6)} 

27T2v^|x|(l  +  e2iKo* sin  ^ / 2TriK0£  sin  ip) 


1  —  erf  I  e  \  2kq£  sin  sin2  ~ 


_ fp*  (y)e^K°(|x|+2^  sin^>) 


2  6' 


7r2\/2^|x|(l  +  e2i/c°^sin^/27rmo^sin'0) 


^21 


I  k0£  sin  -0  sin  2 


/c0  x  oo. 


7T 


(23) 


In  these  asymptotic  expressions  the  principal  contribution  of  the  dimensionless 
quantity  R0£  in,  say,  the  denominator  of  the  integrand  of  (21),  has  the  value  kq£ sin  ip. 
This  implies  that,  for  the  acoustic  far  field,  the  expansion  (20)  is  formally  convergent 
in  the  upper  kq£  sin  ^>-plane  in  the  region  indicated  in  Figure  2,  outside  the  curve 
2tt\k0\£  sinip  =  exp  (— 21m  k0£  sinip).  In  particular,  the  expansion  ceases  to  be  valid  in 
radiation  directions  ip  ~  0,  7r  parallel  to  the  airfoil  span. 
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1.4.  TRAILING  EDGE  NOISE 
1.4.1  Frequency  dependence  of  Green’s  function 

At  each  step  in  the  above  sequence  of  diffraction  problems  the  incident  wave  has  effectively 
been  approximated  by  the  first  term  in  its  asymptotic  expansion  in  powers  of  1/ y/K^ismAp. 
Because  of  this,  as  soon  as  the  approximation  goes  beyond  G  «  G\  +  G]^  +  G^e  the 
neglected  terms  are  formally  of  the  same  order  as  terms  retained  in  G1E,  G% E,  n  >  2.  The 
procedure  we  have  used  corresponds  to  the  retention  of  a  certain  ‘principal  subseries’  of 
terms  from  a  complete  expansion  of  GLE  and  GTe  in  powers  of  1/ y/n0£  sin  ip]  each  member 
of  the  subseries  describes  the  leading  order  effect  of  edge  diffraction  at  each  stage.  The 
procedure  is  familiar  in  quantum  mechanical  calculations  of  multiple  scattering  [32  -  34], 
where  ‘Feynman  diagrams’  are  used  to  extract  convergent  subseries  of  this  kind. 

To  assess  the  accuracy  of  the  various  approximations  the  ratio  |G(x,  y,  tu)/Gi(x,  y,  u>)| 
is  plotted  against  K0£siaip  in  Figure  3  for  typical  cases  6  =  90°,  135°  for  the  three 
approximations: 

{Gi(x,y,u)  +GlE(x,y,v) 

G  i(x,  y,  w)  +  Gle(x,  y,  w)  +  G^E(x,  y>  w) 

Gx(x,  y,  oj)  +  Gle(x,  y,  w)  +  GTE(x,  y,  w), 

where  G\  is  the  half-plane  Green’s  function  in  the  first  line  of  equation  (10).  The  first  is  the 
single  scattering  approximation  (14),  the  second  includes  in  addition  one  scattering  term 
from  the  trailing  edge  (obtained  by  replacing  the  factor  (l  +  e2tK°esm^ /2mK0£sin,ip)  by  unity 
in  the  denominators  of  (22)  and  (23)),  the  third  is  the  multiple  scattering  approximation 
given  by  (22)  and  (23).  The  comparison  is  made  in  the  acoustic  far  field  where 

|x'|  «  |x|  +  i  sin  'ip  cos  9. 

The  function  T(x)  —  g(x)+if(x )  is  conveniently  computed  using  the  rational  approximations 
given  in  §7.3  of  [31],  which  are  quoted  here  for  ease  of  reference: 

r/\_  1  +  0.926a:  .  .  _  1  _ 

*  ^  ~  2  +  1.792®  +  3.104a;2’  9^X)  ~  2  +  4.142a;  +  3.492a;2  +  6.670a;3  ’ 

where  the  absolute  errors  for  0  <  x  <  oo  are  no  larger  than  0.002. 

All  three  approximations  agree  at  high  frequencies  (Ac0^sin^>  >  7).  The  effective 
agreement  between  the  double  and  multiple  scattering  approximations  when  K0isinif) 
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exceeds  about  0.7  suggests  that  the  multiple  scattering,  ‘principal  subseries’  approximation 
is  likely  to  provide  an  excellent  overall  approximation  at  these  frequencies. 

At  very  low  frequencies  the  airfoil  chord  becomes  acoustically  compact,  and  trailing  edge 
noise  will  be  governed  by  the  approximation  (6)  for  Gi(x,y,u>).  The  broken  curve  in  Figure 
4  is  a  plot  of  the  ratio 

((?  i  (x,  y,u)  )compact 
(G1(x,y,w))half_plane 

for  d  =  90°  in  the  acoustic  far  field.  The  solid  curve  is  a  possible  smooth  interpolation  over 
all  frequencies  which  agrees  with  the  high  frequency  multiple  scattering  solution  (dotted) 
for  Koismip  >  0.6.  Similar  interpolations  valid  over  0  <  kq£  <  oo  are  readily  constructed 
for  other  radiation  directions. 

The  finite  chord  of  the  airfoil  can  modify  significantly  the  directivity  of  the  edge  noise 
relative  to  the  familiar  cardioid  pattern  associated  with  a  semi-infinite  airfoil.  This  is 
illustrated  in  Figure  5  for  ip  =  90°  for  k0£  =1,  5,  10  and  50.  The  heavy  curves  are  polar 
plots  of  the  linear  field  strength  (normalized  to  the  same  maximum  value);  in  each  case 
the  light  curve  corresponds  to  the  sin  directionality  for  the  half-plane  (k0£  =  oo).  At 
the  lowest  frequency  shown  in  Figure  5  (k0£  =  1)  the  radiation  peaks  in  directions  normal 
to  the  airfoil,  and  the  directivity  coincides  with  that  for  a  dipole  source  orientated  in  the 
^-direction.  Multiple  lobes  develop  as  the  frequency  increases,  and  the  directivity  tends  in 
an  oscillatory  manner  towards  that  for  the  half-plane.  In  all  cases  there  are  radiation  nulls 
in  the  directions  0  =  0°  and  180°,  respectively  downstream  and  upstream  of  the  airfoil. 

1.4.2  Frequency  spectrum  of  trailing  edge  noise 

The  influence  of  multiple  scattering  on  the  trailing  edge  noise  frequency  spectrum  is 
evidently  independent  of  the  trailing  edge  geometry.  The  prediction  of  the  edge  noise 
spectrum  for  a  finite  chord  airfoil  is  obtained  from  that  predicted  for  an  airfoil  of  semi-infinite 
chord  merely  by  multiplication  by  the  ratio 

G(x,y,a>)  2 

(Gi(x,y,a;))half_plane 

where  G(x,y,u)  is  an  appropriate  interpolation  between  the  low  and  high  frequency 
representations,  such  as  that  depicted  in  Figure  4. 

To  illustrate  the  situation  at  very  low  Mach  numbers,  let  ^(x,  u>)  denote  the  frequency 
pressure  spectrum  of  the  sound  at  the  far  field  point  x  predicted  for  turbulent  flow  from 
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the  edge  of  a  rigid  half-plane.  This  particular  case  has  been  treated  by  many  authors  [3  -  5, 
7  -  11]  by  means  of  the  representation  in  terms  of  the  boundary  layer  wall  pressure  given 
by  equations  (3)  and  (4)  of  §2.  It  will  therefore  be  sufficient  to  present  the  following  result 
from  reference  [14]  for  the  spectrum  of  the  far  field  sound  (|x|  oo)  produced  by  turbulent 
flow  over  one  side  of  the  airfoil 

$oo(x,a;)  «  a0M  sin2(0/2)  sin^^ry^,  M  =  U/c0  «  1,  a0  =  0.035.  (24) 


In  this  formula  5*  and  $pp(oj)  are  respectively  the  boundary  layer  displacement  thickness 
and  the  frequency  spectrum  of  the  wall  pressure  fluctuations  just  upstream  of  the  edge,  and 
L  is  the  span  of  the  airfoil  wetted  by  the  turbulence.  We  shall  use  the  interpolation  formula 
[14]  (based  on  data  collated  by  Chase  [35]) 


(U/6*)%p{uj) 

M)2 


(a 

[M*/U)2  +  a2]l’ 


ap  =  0.12, 


(25) 


where  v *  is  the  boundary  layer  friction  velocity  at  the  edge.  The  broken-line  curve  in 
Figure  6  is  the  corresponding  prediction  of  the  edge  noise  acoustic  pressure  spectrum 
(C//5*)$00(x,cu)/[a0(p0u2)2M(5*L/|x|2)  sin2(0/2)  sinV>]  for  9  =  ^  =  90°. 


The  following  parameter  values  are  typical  of  model  scale  tests  conducted  in  air  and  in 
associated  numerical  simulations  [8  -  11,  24  -  26] 


£_ 

St 


78, 


u£ 

M  =  0.09,  so  that  k0£  =  — 

co 


,,  £  uS* 
M-— -  « 


<5*  U 


7u>5t 

IT 


(26) 


Using  this  to  calculate  k,0£  in  terms  of  u>6*/U,  the  acoustic  pressure  frequency  spectrum  for 
an  airfoil  of  chord  £  is  given  by 


$(x,w)  =  $oo(x,uO 


G(x,y,oQ 

(Gi(x,y,w))half_pl„e 


(27) 


The  normalized  version  of  this  spectrum  is  plotted  as  the  solid  curves  in  Figure  6  for  the  two 
radiation  directions  V  =  90°  and  6  —  90°,  135°.  Significant  departures  from  the  half-plane 
prediction  ^(xjw)  occur  for  oj5JU  <  0.4;  in  particular  spectral  levels  are  reduced  at 
frequencies  below  that  of  the  spectral  peak,  which  is  shifter  to  higher  frequencies. 
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1.5.  CONCLUSION 

Numerical  simulations  and  model  scale  experimental  studies  of  trailing  edge  noise  at  low 
Mach  numbers  frequently  involve  acoustic  frequencies  that  are  sufficiently  small  that  the 
usual  assumption  that  the  airfoil  chord  may  be  regarded  as  semi-infinite  becomes  invalid. 
The  results  presented  in  this  chapter  indicate  that  errors  in  predictions  of  the  edge  noise 
acoustic  pressure  spectrum  based  on  this  model  are  possible  at  lower  and  intermediate 
frequencies,  and  in  particular  in  the  neighborhood  of  the  spectral  peak,  which  tends  to  be 
shifted  to  higher  frequencies  because  of  the  reduced  efficiency  of  edge  noise  production  as 
the  airfoil  chord  decreases  relative  to  the  acoustic  wavelength.  For  the  model  scale  tests 
reported  in  [8,  9,  11]  and  the  numerical  simulations  in  [24  -  26],  the  errors  are  likely  to 
be  significant  when  u5*/TJ  is  smaller  than  about  0.4.  Considerable  departures  from  the 
half-plane  radiation  directivity  occur  at  all  frequencies,  but  especially  at  reduced  frequencies 
K0i  <  10. 

The  approximation  given  in  this  chapter  for  the  edge-source  Green’s  function  may  be 
used  to  solve  a  variety  of  low  Mach  number  aeroacoustic  problems  using  the  general  solution 
given  by  equation  (1).  Of  particular  interest  are  those  involving  the  interaction  of  discrete 
vortex  structures  with  the  leading  or  trailing  edge.  In  general  a  preliminary  calculation 
of  the  incompressible  flow  in  the  immediate  vicinity  of  the  airfoil  must  be  performed  to 
determine  the  time  dependent  distribution  of  the  vorticity.  For  an  impinging  gust,  for 
example,  it  would  be  necessary  to  take  proper  account  of  additional  vorticity  generated  at 
the  surface  of  the  airfoil  (e.g.,  at  the  trailing  edge),  and  also  of  the  evolution  of  the  gust 
vorticity  as  it  convects  past  the  airfoil,  especially  in  the  non-uniform  regions  of  the  mean 
flow  near  the  leading  and  trailing  edges. 
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Figure  2.  The  series  (20)  converges  in  the  upper  complex  n0£ sin  ^-plane  outside  the 
cross-hatched  region. 
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Figure  3.  Comparison  of  single  (•••),  double  (—  —  — )  and  multiple  ( - )  scattering 

approximations  at  0  =  90°,  135°.. 
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Figure  4.  Interpolation  ( - )  between  the  compact  approximation  (second  of  equa¬ 
tion  (6)  with  2/3  =  0)  ( - )  and  the  high  frequency  multiple  scattering 

approximation  (9),  (22),  (23)  (•••). 
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Figure  6.  Comparison  of  the  predicted  edge  noise  spectrum  (24)  for  the  rigid  half¬ 
plane  (U/5^00('x.,u)/[a0(p0vl)2M(6^L/\x\2)  sin2(0/2)  sin  ip]  with  that 
for  an  airfoil  of  finite  chord  £  given  by  (27),  when  ip  =  90°  for  6  =  90°  and 
135°  and  conditions  (26). 
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CHAPTER  2 

ON  THE  HYDROACOUSTICS  OF  A 
TRAILING  EDGE  WITH  A  DETACHED  FLAP 


SUMMARY 

An  analysis  is  made  of  the  production  of  sound  by  low  Mach  number  turbulent  flow 
over  the  trailing  edge  of  an  airfoil  with  a  single  detached  flap.  An  objective  is  to  quantify 
the  separate  contributions  to  the  production  of  sound  by  turbulence  interacting  with  the 
trailing  edge  of  the  airfoil,  the  trailing  edge  of  the  flap,  and  with  the  leading  edge  of  the 
flap.  The  side-edge  noise  of  part-span  flaps  is  not  discussed. 

An  aeroacoustic  Green’s  function  is  derived  for  an  airfoil  of  large  chord  with  a  detached 
flap  at  relative  angle  of  attack  a  (a2  <C  1)  when  the  chord  of  the  flap  is  acoustically 
compact.  Formulae  are  given  for  calculating  the  ‘self-noise’  produced  at  trailing  edges  by 
boundary  layer  instability;  the  efficiency  of  sound  generation  at  the  edge  of  the  airfoil  is 
shown  to  be  typically  at  least  7dB  larger  than  that  produced  at  the  trailing  edge  of  the 
flap.  The  impingement  noise  generated  by  turbulence  interacting  with  the  flap  leading  edge 
is  expressed  in  terms  of  an  equivalent  dipole  source  equal  to  the  fluctuating  flap-lift-force, 
acting  at  a  distance  ip  to  the  rear  of  the  main  airfoil;  dp  is  determined  as  a  function  of 
the  flap  dimensions,  and  does  not  normally  exceed  about  twice  the  width  h  of  the  slot 
separating  the  airfoil  and  flap.  The  proximity  of  the  dipole  to  the  edge  of  the  airfoil  increases 
the  efficiency  of  sound  production  by  a  factor  proportional  to  h/(ipM )  where  M  <C  1  is 
the  characteristic  edge  flow  Mach  number,  and  modifies  the  directivity  of  the  sound.  The 
Green’s  function  can  be  used  with  data  derived  from  direct  numerical  simulations  of  the 
unsteady  hydrodynamic  flow,  and  provides  an  effective  means  of  calculating  the  radiation 
from  a  knowledge  of  the  incompressible  component  of  the  flow  in  the  edge  region. 
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2.1.  INTRODUCTION 

The  efficiency  with  which  sound  is  generated  by  very  low  Mach  number  ‘hydroacoustic’ 
flows  is  so  small  that  the  unsteady  motion  in  the  source  region  is  effectively  indistinguishable 
from  that  of  an  incompressible  fluid.  Great  caution  must  therefore  be  exercised  in 
interpreting  numerical  predictions  of  the  minute  acoustic  byproduct  of  the  flow,  and  it  is 
unwise  to  invest  too  much  confidence  in  direct  numerical  simulations  of  sound  production, 
which  are  often  dominated  by  numerical  rather  than  acoustic  noise.  It  is  particularly 
difficult,  for  example,  to  make  reliable  estimates  of  the  sound  produced  by  turbulence 
in  the  neighborhood  of  a  solid  surface  whose  radius  of  curvature  is  large  relative  to  the 
acoustic  wavelength.  Both  the  amplitude  and  directionality  of  the  sound  are  crucially 
dependent  on  the  relative  phasing  of  very  small  source  components  distributed  over  regions 
spanned  by  many  acoustic  wavelengths;  such  phase  variations  are  largely  indeterminate  by 
conventional  numerical  procedures  when  the  flow  is  effectively  incompressible,  and  spurious 
noise  predictions  can  result  from  the  indiscriminate  use  of  computed  data  in,  say,  Curie’s 
[1]  or  Kirchhoff’s  [2]  surface  integral  representation  of  the  sound. 

This  chapter  is  concerned  with  the  important  special  case  of  the  generation  of  sound 
by  low  Mach  number  turbulent  flow  in  the  vicinity  of  the  trailing  edge  of  an  airfoil  whose 
chord  is  generally  large  relative  to  the  acoustic  wavelength.  Numerical  predictions  in  which 
Curie’s  equation  [1]  (or  the  Ffowcs  Williams  -  Hawkings  equation  [3])  is  used  with  airfoil 
surface  data  derived  from  an  incompressible  numerical  simulation  of  the  flow  are  known 
to  be  incorrect  [4,  5].  Theoretical  methods  for  studying  this  problem  are  usually  based  on 
‘diffraction  theory’  [2,  6-8],  according  to  which  sound  is  produced  when  the  hydrodynamic 
pressure  field  of  turbulence  is  swept  over  the  edge  in  the  mean  flow  and  diffracted  by  the 
edge,  or  on  Lighthill’s  acoustic  analogy  theory  used  in  conjunction  with  a  compact  Green’s 
function  tailored  to  the  trailing  edge  geometry  [2,  9  -  11].  These  methods  are  equivalent 
at  low  Mach  numbers;  their  advantage  is  that  they  unambiguously  identify  fluid-structure 
interactions  at  the  trailing  edge  as  the  source  of  sound,  as  opposed  to  locations  on  the 
extensive,  near  planar  airfoil  surface  upstream  of  the  edge,  where  a  numerical  simulation  of 
the  predominantly  incompressible  flow  would  supply  no  information  about  source  phasing. 
In  other  words,  as  with  the  original  application  of  Lighthill’s  theory  to  free  field  turbulence 
quadrupole  sources,  both  of  these  theoretical  approaches  provide  estimates  of  the  edge  noise 
from  a  knowledge  of  the  incompressible  characteristics  of  the  flow  near  the  edge. 
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Numerical  predictions  of  trailing  edge  noise  that  are  consistent  with  measurements  have 
been  made  by  Wang  [12,  13]  and  Wang  &  Moin  [14]  by  combining  the  results  of  direct 
numerical  simulations  of  incompressible  trailing  edge  flows  with  an  integral  representation 
of  the  sound  involving  the  compact  Green’s  function.  These  calculations  have  hitherto 
treated  the  acoustic  problem  in  terms  of  the  Green’s  function  for  a  trailing  edge  modeled 
by  a  semi-infinite  rigid  plate  [9],  but  the  method  of  conformal  transformation  permits  more 
general  edge  geometries  to  be  considered,  including  airfoils  of  finite  thickness  [2,  5]  and 
finite  chord  (Chapter  1). 

An  important  special  case  arises  for  a  trailing  edge  with  a  detached  ‘flap’  [15  -  17]  (or 
with  the  deployment  of  a  leading  edge  ‘slat’  [18]).  Noise  is  generated  at  the  side-edges  of 
part-span  flaps  and  also  within  the  spanwise  slot  between  the  airfoil  and  flap.  At  very 
low  Mach  numbers,  typical  of  those  encountered  in  underwater  applications,  a  compact 
Green’s  function  can  be  constructed  for  dealing  with  the  production  of  sound  by  side-edge 
sources  [19,  20];  two-dimensional  analytical  models  of  the  influence  of  the  slot  between  the 
deployed  flap  and  the  airfoil  have  been  examined  by  Howe  [21,  22]  for  simplified  aeroacoustic 
distributions. 

However,  to  derive  reliable  analytical  predictions  of  the  intensities  of  noise  produced  by 
general,  three  dimensional  hydroacoustic  source  distributions  in  the  slot  it  necessary  that 
Green’s  function  takes  explicit  account  of  the  doubly-connected  trailing  edge  geometry.  The 
radiated  sound  can  then  be  calculated  from  measurements  or  from  numerical  simulations 
of  the  unsteady  hydrodynamic  flow  in  the  slot.  In  this  chapter  a  formula  is  developed  for 
the  compact  Green’s  function  for  a  flap  of  large  aspect  ratio  in  the  limit  of  small  mean  flow 
Mach  number,  when  the  chord  of  the  flap  is  small  compared  to  the  acoustic  wavelength. 

In  the  first  instance  this  is  done  by  modeling  the  airfoil  by  a  semi-infinite  rigid  plate 
with  a  detached  flap  in  the  form  of  a  rigid  strip  at  a  small  angle  of  attack  to  the  airfoil. 
The  leading  and  trailing  edges  of  the  slot  and  the  trailing  edge  of  the  flap  are  sources  of 
sound  whose  relative  importance  can  be  estimated  from  the  functional  form  of  the  Green’s 
function.  These  conclusions  can  be  extended  to  an  airfoil  and  flap  of  finite  thickness  and 
rounded  edge  geometry,  provided  the  slot  width  is  large  compared  to  the  thickness  of  the 
flap  near  its  leading  edge.  Similarly,  a  correction  can  be  introduced  in  the  manner  described 
in  Chapter  1  to  take  account  of  the  finite  chord  of  the  airfoil. 

The  general  representation  of  sound  produced  at  low  Mach  numbers  by  fluid-structure 
interactions  is  recalled  in  §2.  The  Green’s  function  for  a  flap  at  small  angle  of  attack  is 
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derived  in  §3,  and  approximations  are  obtained  for  its  behavior  near  the  edges  of  the  flap 
and  airfoil.  The  radiation  efficiencies  of  the  different  edges  are  discussed  (§4),  and  leading 
edge,  flap-impingement  noise  is  identified  as  the  most  important  source  in  the  slot.  It  is 
shown  for  this  case  how  the  noise  can  be  represented  in  terms  of  an  equivalent  dipole  source 
(whose  strength  is  equal  to  the  unsteady  lift  on  the  flap)  placed  at  a  suitable  position  to 
the  rear  of  the  airfoil. 
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2.2.  FORMULATION  OF  THE  PROBLEM 
2.2.1  Trailing  edge  noise 

Consider  the  generation  of  sound  at  low  Mach  numbers  by  turbulent  flow  in  the  vicinity 
of  the  trailing  edge  of  an  airfoil.  It  is  first  assumed  that  the  trailing  edge  region  has  the 
idealized  form  illustrated  in  Figure  1.  The  airfoil  is  at  rest  within  a  mean  stream  at  speed  U 
in  the  positive  (^-direction  of  the  rectangular  coordinate  system  (xi,X2,x3),  and  is  modeled 
by  a  semi-infinite  rigid  plate  occupying  the  half-plane  X\  <0,  x^  =  0;  the  detached  ‘flap’ 
consists  of  a  uniform,  rigid  strip  of  width  t  inclined  at  angle  a  to  the  airfoil.  The  slot 
separating  the  airfoil  and  the  leading  edge  of  the  flap  has  width  h.  The  objective  is  to 
determine  the  influence  of  the  flap  and  slot  on  the  generation  of  aerodynamic  sound,  but  no 
account  will  be  taken  of  the  side-edges  of  part  span  flaps. 

Let  v  and  ft  =  curl  v  respectively  denote  the  velocity  and  vorticity,  and  let  c0  be 
the  speed  of  sound.  When  the  mean  flow  Mach  number  M  =  XJ/c0  <C  1  (so  that 
the  convection  of  sound  by  the  mean  flow  may  be  neglected),  the  acoustic  pressure 
p(x,  t)  =  f^°00p(x,co)e~iuJtdio  in  the  far  field  at  position  x  and  time  t  is  given  formally  by  [2, 

5] 

il—Tj  =  y  |5(x,y,w)-(OAv)(y,w)dsy-  1/ jf  f!(y,u)A^(x,y,w)-ndS(y),  |x| ->  00, 

P°  y  *  (28) 
where  p0  is  the  mean  fluid  density  and  v  is  the  kinematic  viscosity.  The  second  integral  is 
over  the  airfoil  surface  S  (with  unit  normal  n  directed  into  the  fluid),  and  can  be  ignored 
when  the  Reynolds  number  is  very  large.  G(x,y,u)  is  the  time  harmonic  acoustic  Green’s 
function,  which  satisfies 

(V2  +  nl)G  =  5(x  —  y),  k0  =  ~,  (29) 

C0 


and  has  vanishing  normal  derivative  dG(x,y]iv)/dyn  =  0  on  S.  When  M  <  1  the  acoustic 
pressure  is  determined  to  an  excellent  approximation  by  equation  (1)  by  approximating  v 
and  ft  in  the  integrands  by  their  values  for  incompressible  flow. 

The  sound  can  also  be  represented  entirely  as  a  surface  integral  over  S  by  introducing  the 
upwash  velocity  Vi,  defined  by  the  Biot-Savart  formula  [23,  24] 


Vi  (x,  t )  —  curl 


r  fl(y,t)d3y 
vs  47t|x  —  y|  ’ 


(30) 
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where  the  integration  is  confined  to  the  nonlinear  region  of  the  hydrodynamic  flow 
outside  the  viscous  sublayers  on  S  [5].  In  this  case 


P(x,w) 

Po 


j£  ^G(x,y,a;)^(y,w)  -  */n(y,w)  A  ^(x,y,u)j -ndS(y),  |x|  -4  oo,  (31) 


To  use  this  formula  a  preliminary  calculation  would  need  to  be  performed  to  determine  the 
unsteady  vorticity  distribution  S7  close  to  the  airfoil.  The  velocity  vr  is  then  calculated  from 
equation  (3);  it  is  precisely  the  ‘kinematic’  velocity  induced  by  the  vortex  field  fl  when  the 
influence  of  ‘image’  vortices  within  S  (i.e.  the  influence  of  scattering  by  S)  is  ignored. 


Equations  (1)  and  (4)  give  alternative,  but  equivalent  representations  of  the  sound  in 
terms  of  the  same  Green’s  function  G(x,  y,u>)  that  has  vanishing  normal  derivative  on  S. 
The  choice  of  which  formulation  to  use  will  in  practice  depend  on  the  details  of  the  problem 
at  hand. 


2.2.2  Acoustically  compact  flap 

Trailing  edge  interaction  noise  is  determined  by  the  contributions  to  the  integrals  in  (1) 
from  the  neighborhood  of  the  edge  region.  The  wavelength  of  the  sound  ~  S/M,  where 
6  is  the  turbulence  eddy  length  scale  near  the  edge.  The  acoustic  wavelength  (~  1/k0) 
will  therefore  be  much  larger  than  the  flap  chord  i  when  M  <  1,  and  the  unsteady  flow 
structures  that  generate  edge  noise  are  then  very  much  closer  than  an  acoustic  wavelength 
from  the  trailing  edge  region.  This  means  that  the  nondimensional  source  distance 
K0\Jy\  +  y\  <.  1,  and  Green’s  function  may  accordingly  be  expanded  in  powers  of  this 
parameter.  The  leading  terms  in  the  expansion  when  the  observation  point  x  lies  in  the 
acoustic  far  field  are  [2] 

G(x,y,w)  =  Gf0(x,y,w)  +  Gi(x,y,w)  +  ...,  |x  -  j/3i3|  ->  oo,  (32) 


where  i3  is  a  unit  vector  in  the  spanwise  (x3-)  direction,  and 


G0(x,y,u>)  = 


_ez«0|x-T/3i3| 

47t|x  -  2/3 i3 


,  Gi(x,y,w)  = 


-^sin5  y>sin  (§)  ^(y^l*-^! 

7iV27ri|x  -  y3i3| 


(33) 


If  r  =  yjx\  +  x\ ,  the  shortest  distance  of  x  from  the  edge  of  the  airfoil  (the  x3-axis), 
then  t/)  =  sin-1(r/|x|)  is  the  angle  shown  in  Figure  1  between  the  direction  of  x  and 
the  edge,  and  (r,  6)  are  the  polar  coordinates  of  x  in  a  plane  rr3  =  constant,  so  that 
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(. xi,x2 )  =  r(cos0,sin0).  The  function  ip*( y)  =  (p*(yi,y2)  of  the  source  position  y  satisfies 
Laplace’s  equation,  and  has  the  simple  interpretation  as  a  velocity  potential  of  ideal, 
incompressible  flow  in  the  anticlockwise  direction  in  the  figure  around  the  edge,  with  zero 
circulation  about  the  flap.  It  is  normalized  by  the  requirement  that 

<p*(y)  p*(y)  =  Vrism(9'/2),  as  r'-^oo,  where  (yi,y2)  =  r'(cos0', sin 0')>  (34) 


where  ^>*(y)  is  the  velocity  potential  of  flow  around  the  rigid  half-plane  x\  <  0,  £2  =  0  [23, 
24], 

The  zeroth  order  component  Go  in  the  expansion  (5)  has  the  same  structure  as  the 
free  space  acoustic  Green’s  function,  and  when  used  in  (1)  or  (4)  yields  contributions  to 
the  acoustic  pressure  that  are  equivalent  to  the  sound  generated  by  free  field  turbulence 
quadrupoles.  The  main  contribution  to  the  radiation  (whose  amplitude  is  larger  by  a  factor 
~  1/M f  l$>  1  than  that  produced  by  the  quadrupoles)  is  generated  by  ‘edge  scattering’ 
and  is  governed  in  a  first  approximation  by  the  second  term  G\  in  equation  (5).  In  the 
high  Reynolds  number  limit,  when  surface  shear  stresses  are  small  relative  to  the  normal 
surface  stresses,  equations  (1)  and  (4)  therefore  yield  the  following  alternative,  leading  order 
representations 


p(x,L0) 


— p0Ko  sins  ipsm(6/2)  f  dip*(y) 


7 r 


\Z2tH\i 


/ 


dy 


(OAv)(y,u;)eiKolx-ij/3i3ld3y 


p0w-\/wcoSin2  i/>sin(0/2) 
7r\/27r|x| 


(35) 


£^*(y)^in(y,^)emo|x  ty3i3|dS(y),  |x|-hx>, 


where  v\n  =  Vi  ■  n  and  d\\/dt  in  (4)  has  been  replaced  by  —iwv  1. 
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2.3.  GREEN’S  FUNCTION 

2.3.1  Calculation  of  <p * 

The  function  <p*(y)  =  <P*(yu  2/2)  is  equivalent  to  the  velocity  potential  of  two-dimensional, 
incompressible  flow  around  the  trailing  edge  region  of  the  airfoil,  and  may  be  determined 
by  the  method  of  conformal  transformation  [23,  24].  In  the  thin  plate  approximation,  the 
airfoil  occupies  the  negative  real  axis  in  the  plane  of  2  =  y\  +  iy^]  the  flap  lies  along  the  ray 
arg  z  =  -a  between  2  =  he~ia  and  z  =  He~ia,  where  H  =  i+h  (see  Figure  2a). 

The  transformation 


maps  the  fluid  region  bounded  by  the  airfoil  and  flap  onto  the  right  half  of  the  (-plane, 
with  the  ‘upper’  and  ‘lower’  sides  (AooO  and  B00O  respectively)  of  the  airfoil  mapping  onto 
the  imaginary  axis,  as  indicated  in  Figure  2b.  The  flap  maps  into  the  interval  (  =  e  16  to 
(  =  e-i£/\/™  of  the  ray  arg  (  =  -e,  where  e  =  a/2  and  m  =  h/H. 

Irrotational  flow  in  the  anticlockwise  sense  around  the  trailing  edge  region  in  the  2-plane 
corresponds  to  flow  past  the  flap  in  the  (-plane  from  (  =  -ioo  to  (  =  -Moo  that  becomes 
parallel  to  the  imaginary  axis  as  |(|  — >  00.  We  can  write 

ip*(  y)  =  Rew((),  (37) 

where  io(()  is  the  complex  potential  of  the  flow  in  the  (-plane.  The  potential  can  be 
calculated  in  closed  form  when  e  =  a/2  =  0;  the  solution  for  e  ^  0  is  found  by  iteration,  by 
writing 

w(C)  =  wo(()  +  «i>i(()  +  •••,  i.e.  <p*(y)  =  y?$(y)  +  e^i(y)  +  •••■  (38) 

The  zeroth  order  problem  for  w0{()  is  illustrated  in  Figure  2c.  The  flap  lies  along  the  real 
axis  between  (  =  1,  1  / s/m\  the  imaginary  axis  forms  a  rigid  barrier  that  may  be  removed 
by  introducing  an  image  flap  on  the  negative  real  axis.  At  large  distances  from  the  flap  in 
the  2-plane  equation  (7)  requires  that  w(z)  —iy/z;  in  the  (-plane  the  flow  is  uniform  and 
parallel  to  the  imaginary  axis  when  |(|  is  large,  and  wq(()  — »•  ^(0  ~  —iVh(. 

The  functional  form  of  u>o(()  with  this  limiting  behavior  can  be  found  by  Sedov’s  method 
([25],  Chapter  2)  to  be  given  by 

».(o  =  -wb£  -  yr=p)  7=f  •  e = co”stant-  <39) 
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The  value  of  p  is  determined  by  the  condition  that  the  circulation  about  the  flap  should 
vanish,  which  yields 

P  =  1  -  ~Sm,l  >0,  ml  =  1  -  m,  (40) 


K(m') 

where  K(m'),  E{rm!)  are  respectively  complete  elliptic  integrals  of  the  first  and  second 
kinds  [26]. 

Higher  order  approximations  to  w(()  can  be  determined  successively  by  imposing  the 
condition  n  •  Vp*  =  0  on  the  flap.  If  £  =  £  +  it],  then  correct  to  0(e)  we  must  satisfy 

d<Pi  d(Po 


for  l<£<l/y/m,  77  =  0. 


dr]  d£  ’ 

It  follows  from  the  symmetry  of  the  problem  that  dpl/d^  is  an  odd  function  of  77,  and 
therefore,  from  (14),  that  p\  represents  the  field  of  a  monopole  source  distributed  over 
1  <  £  <  1/y^m,  77  =  0.  Accordingly,  application  of  the  method  of  images  supplies  the 
solution 

—'/H  /■lA/ro,  2  f  ^  (3  \  dt 


(41) 


«h(C) 


7T 


^I/V,™ln((2  -  C2)  (Vl  -  mi 2 


(42) 


\fl  —  mt2  )  \/t2  —  1  ’ 

which  has  zero  circulation  about  the  flap. 

A  more  convenient  representation  is  obtained  by  considering  the  corresponding  expression 
for  dwi/d(,  obtained  by  differentiation  of  (15).  The  integral  can  be  transformed  into  an 
integral  around  a  closed  contour  in  the  £-plane  that  just  encloses  the  segment  (1,  l/y/m) 
of  the  real  axis.  Cauchy’s  theorem  permits  the  integral  to  be  replaced  by  one  along  the 
imaginary  axis  together  with  a  contribution  from  a  simple  pole  at  t  =  (  (in  Re  £  >  0),  i.e. 

2(y/ff  r'/yM  (1  -  p  -  mt2) 


dwi 

~d( 


7T 


I! 


dt 

\/t2  —  1\/1  —  mt2  (^2  _  C2) 


y/H{l  -  P  -  m(2)  £ VH  [°°  (1-p  +  mX2) 


LI 


dX 


a/1  -  C2\/l  -  ra£2  tt  J- 00  \/l  +  A2\/l  +  mX?  (A2  +  £ 2 ) 
Hence,  using  this  and  equation  (12),  we  find 


(43) 


dw  _  dip*  .dp* 
dz  ~  dyi  dy2 


dw0  dw  1 

+  e- 


dz  ~  dz 

i{z  -  (1  -  mi  /,  .  ia\  a  f~  P  +  (1  -  0)R\  i p 


2 \fz\fz  —  h\fz 


IEL  (i .  _  SL  r 

-H\2J  47 r  /- 


00  \/H  +  p2\/h  +  ]JL2  (fl2  +  z) 


(44) 
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The  real  and  imaginary  parts  of  this  formula  determine  the  functional  form  of  dtp*  (y)  /  dy , 
which  may  be  used  in  the  first  of  equations  (8)  to  evaluate  the  edge  noise  when  the  vorticity 
and  velocity  distributions,  Cl  and  v,  are  known  near  the  flap  (c.f.  [13,  14]). 

To  calculate  the  sound  using  the  second  of  equations  (8),  note  first  that  for  a  thin  airfoil 
it  can  be  written  in  the  form 


p(x,w) 


poujy/m^ sin ^  V>sin(fl/2) 
7r-\/27r|x| 


f  [»>,(y)]«i„(y,w)e“"l’I-iw'*l«lS(y))  |x|  ->  oo,  (45) 

J  s+ 


where  the  integration  is  over  the  ‘upper’  surfaces  S+  of  the  airfoil  and  flap,  and  [</?*]  =  <p*+—p>*_ 
is  the  jump  in  the  value  of  <p*  in  passing  from  the  point  y  on  the  lower  surface  to  the 
corresponding  point  on  the  upper  surface.  We  can  determine  [</?*]  by  integrating  the  jump 
d[w\/dz  defined  by  equation  (17).  Introducing  the  approximation  1  +  ia/2  «  e'f,  we  find 

[p*]  =  2^/lyi^~yi)  ~  2H^£  (tan-1  Jh)  X) 


-(1  -  0) T  ^tan-1  (\J\yi\/h^j  ,  m'^j  j,  yx  <  0  on  the  airfoil, 


(46) 


where 


[</?*]  =  2^  — — — —  -  2 S^E  ^sin  1  (^H{y\\  -  W/v^)  > m' ) 

-(1  -  0)F  ^sin-1  (^JW{y\\  -  h)/y\\Zj  ,  m’^j  |,  h  <  y\\  <  H  on  the  flap,  (47) 


£{£,m')=  [  Jl-  rri'  sin2  /i  d/i,  T(£,m')  =  [  ,  = 

Jo  Jo 


(48) 


—  m'  sin  [jl 

and  y\\  denotes  distance  measured  parallel  to  the  surface  of  the  flap  (i.e.,  in  a  direction 
inclined  at  the  flap  angle  a)  from  the  edge  O  of  the  airfoil. 


2.3.2  Critical  behavior  near  the  edges  of  the  flap  and  airfoil 

The  solid  curves  in  Figure  3  represent  the  potential  flow  streamlines  V’o  =  constant  for 
the  zeroth  approximation  wo  =  <Pq  +  q  when  i  =  Ah  (m  =  0.2),  calculated  from  equation 
(12)  by  setting  £  =  \fzh.  The  volume  flux  between  neighboring  streamlines  is  the  same 
in  all  cases;  thus  convergence  of  the  streamlines  near  the  leading  and  trailing  edges  of 
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the  slot  between  the  airfoil  and  flap,  and  near  the  trailing  edge  of  the  flap  correspond  to 
regions  where  is  large  and  where,  according  to  equation  (8),  the  efficiency  of  edge 
scattering  of  turbulence  sources  becomes  large.  The  stagnation  streamline  intersects  the 
flap  symmetrically  and  corresponds  to  the  broken  line  curve  in  the  figure. 

The  behaviors  of  <p*( y)  =  Re  w(()  near  the  edges  of  the  flap  and  airfoil  are  determined 

by  expanding  w(()  =  w{Jzjh)  about  the  critical  points  at  z  =  0,  h,  H  (=  t  -f  h).  To 

the  present  order-e  approximation,  the  singularities  at  2  =  0,  h,  and  H  occur  respectively 
at  the  edge  of  the  airfoil,  and  at  the  leading  and  trailing  edges  of  the  flap.  By  expanding 
the  second  line  of  (17)  about  each  of  these  points  and  integrating  with  respect  to  z  we  find 
(after  discarding  irrelevant  constants  of  integration) 

airfoil  trailing  edge:  w  ~  —  i  _  flyfi  +  | 

flap  leading  edge:  w  ~  \fz  —  heT  [(1  —  /3)  (l  +  f ) 

flap  trailing  edge:  w  ~  —  iyj~z  —  (h  4-  £)e^ / 3\Jl  +  j 

The  corresponding  approximations  for  <p*(y )  —  Re  w  can  be  written 

airfoil  trailing  edge:  <p*  ~  y/ro  sin  (1  —  +  \ 

flap  leading  edge:  <p*  ~  ^cos  (^)  [(1  -  /?)  (l  +  f)  -  l]  Q 

flap  trailing  edge:  ip*  ~  ^/FffSin  /?^/l  +  f 


where  (ro,#o)>  (rfc,0/ i)>  (rH,  6h)  are  respectively  the  polar  coordinates  of  the  source  point 
(z/1,2/2)  relative  to  the  edge  O  of  the  airfoil,  to  the  flap  leading  edge,  and  the  flap  trailing 
edge. 
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2.4.  TRAILING  EDGE  AND  FLAP  GENERATED  SOUND 


2.4.1  Trailing  edge  noise 

An  important  component  of  the  sound  generated  by  the  trailing  edge  region  can  be 
attributed  to  turbulence  whose  characteristic  length  scale  is  small  compared  to  the  flap 
chord  t.  In  this  case  the  effective  acoustic  sources  at  the  edges  labelled  0,  A  and  B  in 
Figure  4a  will  be  statistically  independent,  and  their  relative  contributions  to  the  far  field 
sound  can  be  estimated  by  inserting  the  corresponding  approximations  (23)  for  <p*(y)  into 
the  second  of  equations  (8)  and  performing  local  surface  integrations. 

The  upwash  velocity  vt  is  produced  by  the  turbulence  convecting  past  the  edges, 
augmented  by  the  induced  velocity  of  any  additional  vorticity  shed  from  the  trailing  edges 
(to  satisfy  the  unsteady  Kutta  condition  [27]).  The  effect  of  this  shedding  is  greatly  to 
reduce  the  importance  of  the  edges  O  and  B  as  sources  of  sound,  because  the  induced 
velocity  of  the  shed  vorticity  acts  to  oppose  that  produced  by  impinging  turbulence.  This 
does  not  occur  at  the  leading  edge  A:  provided  the  flow  does  not  separate  near  A,  any 
vorticity  generated  there  is  swept  downstream  over  the  surface  of  the  flap,  and  its  induced 
velocity  is  therefore  cancelled  by  an  equal  and  opposite  velocity  produced  by  image  vortices 
in  the  flap.  Thus,  impingement  noise  generated  at  A  is  likely  to  be  the  principal  source  of 
sound  in  the  trailing  edge  region. 


However,  it  is  possible  for  the  trailing  edges  O  and  B  to  remain  significant  sources  of  high 
frequency  sound.  It  follows  from  equation  (8)  and  by  inspection  of  (23)  that  the  relative 
intensities  of  the  corresponding  acoustic  pressures  po,  Pb,  say,  is  given  by 


Po 

,Pb, 


Uo\2  l  (1  -  pf 
vB)  h  P2 


(51) 


where  Vo,  vB  are  characteristic  turbulence  velocities  (e.g.  friction  velocities)  respectively 
near  the  trailing  edges  of  the  airfoil  and  flap.  The  relative  efficiency  factor 
10  x  log10[(£/h)(l  -  P)2/P2]  (dB)  determines  the  ratio  of  the  intensities  of  the  sound 
produced  at  O  and  B  when  the  magnitudes  of  the  turbulence  velocities  are  nominally  equal; 
it  is  plotted  as  a  function  of  £/h  in  Figure  5,  from  which  it  is  seen  that  the  efficiency  of  edge 
noise  production  at  the  trailing  edge  O  of  the  airfoil  exceeds  that  at  the  trailing  edge  of  the 
flap  typically  by  7  dB  or  more. 
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2.4.2  Equivalent  flap-impingement-noise  dipole 

Equations  (8)  and  the  second  of  (23)  show  that  the  acoustic  pressure  generated  by  small 
scale  turbulence  impinging  on  the  leading  edge  of  the  flap  at  A  is  given  by 

fp0u)y/2iK^sin^  ip  sin(6 /2)elK°\x\  r°°  re+h 


p(x,w) 


3|  I 
7T2  X 


[  dy3  j  Jy\\  -  hv±(y\\,y3,u)e  iKoCOS^3dy||, 

J  — oo  Jh  v 

|x|  — >  oo,  (52) 


where  v±(y\\,y3,u)  is  the  normal  component  of  the  upwash  velocity  (directed  ‘upwards’ 
from  the  flap),  expressed  in  terms  of  the  distance  y\\  measured  along  the  flap  from  A,  and 


/  =  / 


(53) 


The  turbulence  impinging  on  the  leading  edge  A  produces  lift  fluctuations  on  the  flap. 
The  magnitude  of  this  lift  force  will  be  influenced  to  some  extent  by  the  proximity  of  the 
airfoil,  but  for  small  scale  turbulence  the  influence  will  be  small,  and  the  force  will  be  close 
to  its  value  when  the  airfoil  is  absent.  The  lift  per  unit  span  at  frequency  u>  for  a  thin- plate 
flap,  when  the  reduced  frequency  u>£/U  >  1,  and  when  the  influence  of  the  main  airfoil  is 
ignored,  is  given  by  [2,  28] 

F(y3,u)  =  -2ip0uj\/i  J  yfy\\  -  ~hv±(y\\,y3,u})dy\\.  (54) 

This  is  equivalent  to  the  classical  Sears  formula  [29]  when  v±(y\\,y3,u>)  is  assumed  not  to 
depend  on  the  spanwise  coordinate  y 3 

Let  us  determine  the  sound  generated  by  this  force  when  it  is  assumed  to  be  distributed 
along  a  spanwise  line  at  distance  Ip  to  the  rear  of  the  airfoil  in  the  absence  of  a  flap 
(see  Figure  6).  The  force  on  the  fluid  is  (0,  —  F,  0)  and  the  acoustic  pressure  is  therefore 
determined  by  [2] 

(V2  +  k l)p  =  -J^-(6(x  1  -  £F)8(x2)F{x3,wj),  (55) 

subject  to  dp/dx 2  =  0  on  the  airfoil.  The  solution  can  be  expressed  in  terms  of  the 
component  (?i(x,  y,o>)  of  Green’s  function  (5),  where  </?*(y)  is  given  by  (7)  in  the  absence 
of  the  flap: 
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r\ 

p(x,w)  =  -  j  G1(x.,y,u)  —  (6(y1  -  £F)S(y2)F(y3,w))di y 


-  sin  2  ip  sin  (0/ 2)  elKo  lx'  f°° 


2ttV2^|x| 


f°°  F(yz,u)e-iK°cos^dy3 

J  —  OO 


I  t  pQcj\/2m0  sin2  ^sin(fl/2)em°|xl 
4^f  7T  2  |x| 

x  f  dy3  f  Jy\\  v1_(yi\,y3,ui)e~lKoCOS'l)y3dy\\, 

J  —  oo  */  h 


x|  — »■  oo.  (56) 


This  will  coincide  with  the  directly  calculated  radiation  (24)  if  lF  —  £/4/2.  In  other 
words,  provided  the  dipole  F  is  situated  downstream  of  the  airfoil  trailing  edge  0  when  the 
presence  of  the  flap  is  ignored  at  a  distance 


If  _  (W 

h  4  [(1  - /?)  (l  +  f )  -  if 


This  function  is  plotted  in  Figure  6.  Evidently  tF/h  increases  very  slowly  as  the  chord  of 
the  flap  increases  (ultimately  proportional  to  [ln(16^//i)/4]2),  and  iF  does  not  exceed  twice 
the  slot  width  h  until  ijh  >  12. 

When  the  the  statistical  properties  of  F(y3,  u)  are  known  the  frequency  spectrum  $(x,w) 
of  the  acoustic  pressure  at  x  (which  satisfies  (p2(x,  t))  —  4>(x,  to)du>,  where  (  )  denotes 

an  ensemble  average)  is  readily  calculated.  For  stationary  random  turbulence 


(F(y3,  uo)F*(y'z,  to1))  =  5  (to  —  to')^FF(to)7Z(\y3  —  y'z\,to),  (58) 


where  $FF(co)  is  the  frequency  spectrum  of  the  lift  per  unit  span  (defined  such  that 
(|_F(y3,u;)|2)  =  /f^4>FF(k>)dw)>  and  F,(y,to)  is  the  corresponding  spanwise  covariance  (for 
which  71(0,  to)  =  1). 


In  terms  of  these  definitions,  it  is  easily  deduced  from  the  second  of  equations  (29)  that 

,2 /„ . . 

|x|  — »•  oo,  (59) 


sin  sin2  (0/2)  L£3  .  . 

4(x-w)  K 


where  £3  —  ff°  7Z(y,u>)dy  is  the  spanwise  correlation  length  of  the  lift  and  L  »  £3  is  the 
span  of  the  flap  wetted  by  the  turbulent  flow.  The  spanwise  correlation  scale  is  frequency 
dependent,  and  typically  £3  ~  F j  to,  so  that  k0£3  ~  M  1;  this  condition  has  been  used  in 
deriving  (32).  The  order  of  magnitude  of  the  flap  mean  square  acoustic  pressure  is 
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which  is  characteristic  of  edge  generated  sound  when  the  airfoil  chord  is  large  compared  to 
the  acoustic  wavelength. 

2.4.3  Airfoil  and  flap  of  finite  thickness 

These  detailed  conclusions  are  applicable  to  the  more  general  trailing  edge  configuration 
sketched  in  Figure  4b.  The  formulae  (23)  for  <p*( y)  now  constitute  ‘outer  approximations’ 
to  the  potential  function  defining  irrotational  flow  around  the  trailing  edge  region.  <p*(y) 
varies  in  a  more  complicated  manner  close  to  the  rounded  edges  of  the  flap,  or  near  the 
thickened  geometry  of  the  airfoil  trailing  edge,  but  will  tend  to  the  forms  given  in  (23)  when 
the  distance  from  an  edge  exceeds  the  edge  radius  of  curvature.  This  indicates  that  the 
above  predictions  of,  say,  the  flap  impingement  noise  should  still  be  valid  for  components  of 
the  upwash  velocity  whose  length  scales  are  larger  than  the  leading  edge  radius  of  curvature 
of  the  flap.  But  these  would  be  expected  to  be  the  dominant  small  scale  noise  sources, 
since  the  efficiency  of  sound  generation  tends  to  fall  off  rapidly  when  the  turbulence  scale 
becomes  smaller  than  the  radius  of  curvature  [2].  In  particular  the  result  in  Figure  6  for  the 
location  of  the  effective  flap  dipole  source  should  remain  true  for  the  dominant  turbulence 
sources. 
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2.5.  CONCLUSION 

The  sound  produced  by  a  trailing  edge  flap  with  a  detached  flap  has  three  principal 
sources.  First  there  is  the  ‘self  noise’  attributed  to  small  scale  turbulence  produced  by 
instabilities  in  the  boundary  layers  upstream  of  the  airfoil  and  flap  trailing  edges.  This 
is  of  relatively  high  frequency,  because  the  sound  that  would  be  generated  by  larger  scale 
turbulence  eddies  (large  on  a  scale  of  airfoil  or  flap  thickness  near  the  edge)  tends  to  be 
cancelled  by  that  produced  by  vorticity  shed  from  the  edge  in  accordance  with  the  Kutta 
condition.  Second,  impingement  noise  is  produced  when  turbulence  is  swept  past  the 
leading  edge  of  the  flap,  creating  an  unsteady  lift  force  that  generates  sound  by  interaction 
with  the  trailing  edge  of  the  airfoil.  Finally,  the  side  edges  of  part-span  flaps  constitute  an 
important  source  of  lower  frequency  noise,  associated  with  the  interaction  of  side-edge  lift 
vortices  with  the  side-edge  and  trailing  edge. 

The  formulae  given  in  this  chapter  for  the  compact  Green’s  function  for  a  trailing  edge 
with  a  single  detached  flap  can  be  used  to  estimate  the  trailing  edge  self-noise  and  the 
flap  impingement  noise  in  flows  at  low  Mach  number.  The  analytical  results  for  a  flap  at 
a  modest  angle  of  attack  a  (where  a2  <C  1)  indicate  that  the  self-noise  produced  at  the 
edge  of  the  airfoil  is  typically  at  least  7dB  in  excess  of  that  produced  at  the  flap  trailing 
edge.  The  strength  of  the  noise  generated  by  impingement  on  the  flap  leading  edge  can 
be  calculated  by  first  determining  the  unsteady  lift  F(u)  per  unit  span  produced  by  the 
turbulence  when  the  flap  is  regarded  as  an  isolated  airfoil,  for  example,  by  means  of  a 
convenient  ‘strip-theory’  approximation  based  on  Sears’  formula  for  the  lift  produced  by  a 
gust.  The  radiation  is  then  equal  to  that  produced  by  a  dipole  of  strength  —  F  situated  at 
distance  ip  to  the  rear  of  the  edge  of  the  main  airfoil;  this  distance  is  given  as  a  function 
of  flap  dimensions  by  equation  (26).  The  dipole  radiation  intensity  is  increased  from  its 
classical  free  field  value  by  the  proximity  of  the  airfoil  trailing  edge,  typically  by  a  factor 
~  0(1  /M)  (^>  1);  similarly  its  directivity  is  modified  by  the  presence  of  the  airfoil. 

In  this  chapter  it  has  been  assumed  that  the  airfoil  chord  is  much  larger  than  all  other 
length  scales,  including  the  acoustic  wavelength.  This  condition  can  be  relaxed,  however 
(provided  the  characteristic  acoustic  wavelength  remains  large  compared  to  the  flap 
chord),  without  changing  the  principal  conclusions,  by  introducing  the  frequency-dependent 
correction  factor  given  in  Chapter  1  for  dealing  with  airfoils  of  finite  chord.  Predictions  of 
this  chapter  can  therefore  be  used  to  validate  acoustic  numerical  prediction  schemes  applied 
to  a  finite  chord  airfoil  with  a  detached  flap  of  suitably  simple  geometry. 
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Figure  2.  (a)  Thin  plate  airfoil  and  flap  in  the  z-plane;  (b)  Image  in  the  £-plane; 
(c)  Configuration  in  the  £-plane  in  the  zeroth  approximation. 
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Figure  4.  (a)  Turbulence  interacting  with  the  airfoil  and  flap;  (b)  Generalized  trailing 
edge  configuration. 
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Figure  6.  Line  of  action  of  equivalent  dipole  lift  force  F  for  calculating  flap  leading- 
edge-impingement  noise. 
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